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ABSTRACT 

V 

c=~A  comprehensive  treatment  of  the  classical  decay  rates  for  a  molecule  in 
the  vicinity  of  a  spherical  surface  is  presented  through  the  application  of  the 
work  of  van  del  Pol  and  Bremmer  and  that  of  Fock.  This  theory  takes  full 
advantage  of  the  Hertz  vector  formalism,  which  is  mathematically  simpler  than  the 
widely-adapted  Lorenz-Mie  approach  which  uses  the  complicated  vector  harmonic 
expansions.  Results  are  obtained  for  both  radiative  and  nonradiative  transfer 
when  the  molecule  is  located  outside  or  inside  the  surface.  Numerical  results 
are  given  for  the  cases  of  a  surface,  cavity  and'^erosol^  and  comparison  with 
other  works  are  made.  .  -  -<• 
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I.  INTRODUCTION 

Ever  since  the  role  of  surface  roughness  was  recognized  as  prominent  in 

leading  to  such  dramatic  surface  phenomena  as  surface-enhanced  Raman  scattering 

and  other  photochemical  processes,^  the  problem  of  the  decay  rates  for  molecules 

2-5 

in  the  vicinity  of  a  spherical  particle  has  been  investigated  intensively. 

While  the  classical  electromagnetic  approach  has  been  folloved  to  establish  the 

2 

results  for  the  reduced  (normalized)  decay  rates  in  both  the  electrostatic  limit 
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and  the  more  exact  electrodynamical  treatment ,  we  feel  that  a  complete 
analysis  of  such  a  problem  has  not  yet  been  available  in  the  literature.  Given 
two  possible  orthogonal  orientations  of  the  molecule  (radial/ tangential),  the  two 
possible  locations  of  the  molecule  (outside/ inside  the  sphere)  and  the  two 
different  kinds  of  energy  transfers  during  the  decay  process 
(radiative/nonradiative),  there  are  altogether  eight  problems  to  be  solved.  In  a 
paper  following  the  analogous  treatment  of  Chance,  Prock  and  Silbey  for  a  flat 
surface,6  Ruppin  solved  the  complete  problem  for  a  molecule  located  outside  the 
sphere.  Very  recently.  Chew4  published  results  using  the  energy  flux  method'*' 
for  both  the  cases  where  the  molecule  is  located  outside  and  inside  the  sphere, 
and  the  equivalence  with  the  results  obtained  from  Green's  dyadic  method6  is 
proven  in  the  limit  where  all  the  dielectric  constants  are  real.  Hence,  all  the 
nonradiative  transfers  are  ignored  in  Chew's  treatawnt.*  As  a  matter  of  fact, 
Eqs.  (6)  and  (7)  in  Ref.  4  are  identical  to  the  results  for  radiative  transfer 
for  a  molecule  located  outside  the  sphere  as  given  is  Eqs.  (27)  and  (28)  in 
Ruppin' s  work.  Hence,  combining  the  work  of  Ruppin  and  Chew,  we  conclude  that 
six  of  the  eight  problems  are  solved,  leaving  the  two  problems  of  nonradiative 
rates  for  a  molecule  inside  the  sphere  to  be  solved.  Moreover,  all  these  past 


treatments  have  been  based  on  the  Lorenz-Mie  theory  established  by  the  expansion 


of  the  field  quantities  using  vector  spherical  harmonics. 
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In  this  present  work,  we  present  a  complete  solution  of  the  eight  problems 

via  a  different  approach,  namely,  the  Herts  vector  formalism. This  formalism 

has  a  mathematically  simpler  feature,  since  the  expansions  of  the  Herts  functions 

are  more  of  a  scalar  type,  and  one  can  see  easily  from  the  Green's  function 

property  that  there  exists  a  simple  transformation  which  relates  the  fields  for 

the  outside-molecule  case  to  those  where  the  molecule  is  located  inside  the 

sphere.  Hence  it  is  not  necessary  to  solve  again  the  boundary  value  problem  for 

4 

the  inside -molecule  case  as  is  done  in  Chew's  work.  We  elaborate  the  theory  in 

Sec.  II,  where  we  shall  see  that  for  the  case  where  the  molecule  is  inside  a 

sphere  with  a  complex  dielectric  constant,  a  straightforward  application  of  the 

3  4 

methods  of  Rupp in  and  Chew  is  not  appropriate,  such  that  one  must  resort  to 
more  microscopic  models  for  the  treatment  of  the  bulk  decay  rates  (i.e.,  in  the 
absence  of  the  surface).  Numerical  examples  are  given  in  Sec.  Ill,  where  we 
shall  illustrate  that  the  neglect  of  the  nonradiative  rates  as  in  Chew's 
formalism  can  lead  to  serious  errors  in  some  cases.  Discussion  and  conclusion  of 
the  present  work  are  given  in  Sec.  IV. 


II .  Theory 

It  is  well  known  that  the  electromagnetic  scattering  problem  in  the 
presence  of  a  dielectric  sphere  csn  in  general  be  formulated  in  two  different 
ways.  While  the  scattering  of  a  plane  wave  from  a  sphere  has  been  well  treated 
in  both  the  Lorens -Hie  approach^  and  the  Herts  vector  (Debye  potential) 

g 

approach,  it  seems  that  the  problem  of  dipole-sphere  interactions  (hence  the 
problem  of  aiolecular  lifetimes)  has  always  been  treated  by  the  former  approach 
using  the  vector  harmonics  expansion. Here  we  want  to  solve  the  lifetime 
problem  in  the  Herts  vector  (ft)  formalism  using  expansions  in  ordinary  spherical 
functions.  For  those  mathematical  results  which  have  been  obtained  already  by 


4 


the  former  approach,  we  shall  only  indicate  briefly  how  the  same  results  may  be 

obtained  in  our  present  approach,  with  possible  slmplif ications  being  emphasized. 

Por  a  dipole  located  outside  the  sphere,  the  problem  has  been  solved  for 

9  10 

both  the  vertical  (radial)  and  the  horizontal  (tangential)  dipoles.  Assuming 
nonmagnetic  media  and  the  region  outside  (inside)  the  sphere  (of  radius  e)  to  be 
characterized  by  the  local  dielectric  function  (£2),  the  Herts  vectors  (along 
the  radial  direction)  can  then  be  expressed  in  the  following  series  expansions  in 
terms  of  the  various  "Mie  coefficients": 


For  a  vertical  dipole  located  at  (d.0.0)  with  ■  k^d,  ■  k^a  and  k^  ■ 


C.W/C, 


r  ikj  £  (2n+l)  h^1)(y1)Ijn(k1r)  +  »nh^1)(k1r) l^cose)  ,  a  <  r  <  d 


ikj  £  (2n+l)ljn(yi)  ♦  »nh^l)(yi)l  hJ^Ujc)  Pn(cose)  ,  d  < 


■  -I  l  TM  W)rn<“*8> 

Ka  n  a 


T  <  ft 


2  n-0 


where  Jn  and  h^  '  are  the  ordinary  spherical  Bessel  and  Henkel  functions.  The 
nonvanishing  field  components  are  then  given  by 


rd  sine  ae  (tln8  ae*  • 


1  »  / 
rd  arae  *rn*  • 
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For  a  horizontal  dipole  located  at  (d.0.0).  wo  need  two  Herts  vectors  with  one 


along  r  and  one  along  fi,  given  by  (i  *1,2) 
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where  P,  Q  are  expressed  in  series  fom  as  (with  (x)  S  xj  (x)  and  C  (x)  » 

n  n  n 

xh(1)(x)) 
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The  nonvanishing  field  components  are  then  given  by 
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The  "Hie  coefficients"  in  the  above  equations  are  defined  as 


*- - — JiT. 
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B  -  c2^n^p2)*n^pl > 

'  clhn1>(ol)^<p2) 


Cn  '  V*  W  ■  ^(“iW 
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From  this  formulation  and  by  making  use  of  the  veil-known  properties  of  Green's 


functions,  ve  can  easily  obtain  via  a  simple  transformation  T  the  solution  for 


the  reflected  and  transmitted  fields  for  the  case  where  the  dipole  is  located 


inside  the  sphere  without  having  to  solve  again  the  boundary  value  problem.  If 


ve  define  T  by 


kj,tj, . . . 


^2*^2*  *  *  * 


^2  *  e2 *  *  *  * 


kj ,tj , « . . 


then  it  is  not  difficult  to  check  that  the  "inside  solutions"  are  easily  obtained 


from  the  oustide  solutions  through  the  operation  with  T. 


Now  let  us  apply  the  above  results  for  our  decay  rate  problem.  According 


to  the  classical  approach,  there  are  two  different  methods  for  calculating  the 


decay  rate  of  the  molecule  as  Induced  by  the  presence  of  the  surface.  In  the 


Green's  dyadic  approach,  the  dipole  is  modeled  as  a  damped  harmonic  oscillator 


being  driven  by  the  surface  field.  The  induced  decay  rate  is  then  given  by 


Y  •  Y0  ♦  -  la  G  , 
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where  Yq  is  the  decay  rate  in  the  absence  of  the  boundary  surface  (i.e.,  with  the 
whole  space  filled  with  the  dielectric  medium  in  which  the  dipole  is  located),  a 
■  e  /m,  uj  being  the  emission  frequency  (which  is  assumed  to  change  negligibly)* 
and  Im  G  is  the  imaginary  part  of  the  dyadic  function  defined  as  the  surface 
field  (Eg)  reflected  at  the  dipole  site  per  unit  dipole  moment  (jl).  The  other 
approach  would  be  to  calculate  the  rate  of  the  energy  carried  away  from  the 
molecule,  which  can  further  be  divided  into  that  radiated  to  infinity  (radiative 


shall  apply  the  fields  obtained  from  Eqs.  (1)-(10)  to  calculate  the  y's  given  in 
Eqs.  ( 13)-( 15)  for  the  following  three  cases: 
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Case  (i):  The  surface  problem 

Here  p  is  located  in  vacuum  (c^  ■  1)  outside  a  sphere  of  dielectric 

X  A  2  3 

constant  €2  ■  £*  +  ie",  and  Yq  is  then  given  by  ^ckjp  /W  *  2ak1/3w.  Using  Eqs. 

4 

(1),  (3a)  and  (13),  we  obtain  for  a  vertical  dipole 


1  r  n  *  y  1  *  2 

1  +  |  Re  £  ( 2n+l )n( n+l)Bn[ 


Furthermore,  using  Eqs.  (l)-(3),  the  Poynting  vector  can  readily  be  calculated, 

3  4 

and  from  Eq.  (15)  we  obtain  ’ 
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Similarly,  using  Eqs.  (4)-(10),  (13)  and  (15),  we  obtain  the  corresponding 

results  for  a  horizontal  dipole 
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(20) 


The  nonradiative  rates  can  then  be  obtained  from  the  difference  between  the  total 
and  radiative  rates. 

Case  (ii)s  The  cavity  problem 

Here  p  is  located  inside  a  spherical  vacuum  (cj  *  1)  with  a  dielectric  (e^ 
■  e'  +  ie")  occupying  the  whole  space  outside  the  vacuum.  The  total  decay  rates 
can  readily  be  obtained  from  the  outside  solutions  (Eqs.  (17)  and  (19))  through 
the  transformation  T  (Eq.  (12))  as  discussed  above,  whereby  we  obtain** 


XL 
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(22) 


with  E  and  F  given  by 
n  n 

E  -  T{B  }  , 

n  n 

(23) 

F  -  T{A  }  . 

n  n 

(24) 

He  want  to  remark  that  in  this  configuration  the  decay  rates  can  only  be  purely 
radiative  or  purely  nonradiative,  depending  on  whether  the  outside  medium  is 
transparent  (c"  -  0)  or  dissipative  (c"  *  0)  with  respect  to  the  molecular 
emission  frequency.  Furthermore,  it  is  interesting  to  note  that  in  the 
perfectly-conducting  limit  ( | £2 | » | | )  where  we  have* 
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both  Eqs.  (21)  and  (22)  vanish.  Thus,  in  the  present  classical  treatment,  the 


molecular  lifetime  becomes  infinite  in  a  perfectly-conducting  cavity.  This  is 


understandable  since  in  this  approach  only  two  mechanisms  can  "cause"  the 


molecule  to  decay,  i.e.,  either  for  energy  being  brought  to  infinity  (radiative 


transfer)  or  dissipated  into  a  host  medium  (nonradiative  transfer).  Since  no 


field  can  penetrate  into  the  perfect  conducting  environment  and  hence  both  these 


transfers  cannot  occur,  it  is  not  surprising  to  see  that  both  decay  rates  vanish 


in  Eqs.  (21)  and  (22).  In  the  more  exact  quantum  treatment,  however,  the 


vanishing  of  the  decay  rates  in  a  perfectly-conducting  cavity  can  only 
. .  12 


restrictively  occur. 


Case  (iii)s  The  "aerosol"  problem1 


Here  m  is  located  inside  a  dielectric  sphere  (C2  *  e’  +  ie")  with  the 


outside  being  vacuum  (c^  ■  1).  He  shall  further  divide  this  into  the  following 


two  cases  for  discussion: 


(e)  Transparent  dielectric 


This  case  has  already  been  treated  by  Chew.  In  the  present  approach,  the 


total  decay  rate  can  again  be  obtained  directly  from  Eqs.  (17)  and  (19)  through 
T.  With  yn  being  given  by  ^ckjcj^W^  ■  2ak^/3ue9,*  we  obtain  the  same 


expressions  as  given  in  Eqs.  (21)-(24).  The  radiative  transfers  in  this  case  can 
similarly  be  obtained  as  in  Eqs.  (18)  and  (20)  via  the  transformation  ■ 


niout)  to  give* 
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(26) 


(27) 


In  this  case,  since  is  real  (transparent,  e"  -  0),  there  is  no 
dissipative  (nonradiative)  effects,  and  it  has  been  shown  explicitly  by  Chew  that 
the  radiative  rates  in  Eqs.  (26)  and  (27)  are  identical  to  the  total  rates  as 
given  in  Eqs.  (21)  and  (22). 


(b)  Dissipative  dielectric 

This  is  the  case  which  has  not  been  discussed  before  in  the  context  of 

"surface  problems".  Nevertheless,  the  same  problem  for  a  molecule  in  the  bulk  of 

14 

an  adsorptive  medium  has  received  considerable  attention  in  the  literature. 
The  difficulty  lies  in  the  evaluation  of  Yq  in  Eq.  (13).  Obviously,  the  decay  in 
a  dissipative  medium  of  infinite  extent  can  only  be  nonradiative.  However,  If 
one  applies  Eq.  (16)  directly  to  calculate  Yq  in  this  case,  one  will  get 
divergences  in  the  integral  due  to  the  predominance  of  the  near  field  in  this 
case,  unless  one  adopts  a  cutoff  volume  (forbidden  volume**)  in  the  lower  bound 
of  the  integral.1®  Moreover,  this  approach  has  been  criticised  by  Agranovich  and 
Dubovskli  (AD),1®  who  proposed  a  more  exact  microscopic  treatment  which  shows 
that  the  result  from  such  macroscopic  theory1®  can  be  correct  only  under  certain 
restrictive  conditions.  In  the  following,  we  shall  make  use  of  the  results  of  AD 
to  investigate  the  surface  effects  in  Eq.  (13). 
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According  to  AD,  under  the  condition  that  the  "Mossbauer-type"  effect  (in 
the  optical  region)^  can  be  neglected,  Yq  can  be  expressed  in  a  form  which  is 
very  similar  to  that  obtained  from  Feinberg's  theory,^ 


—  -  -&r  e'V3 

Y  e.!  A  0 
sp  64ir 


where  YSp  is  the  spontaneous  decay  of  the  free  molecule,  Nq  is  the  molecular 
number  density  of  the  dielectric  medium,  A  is  the  emission  wavelength,  and  Z  is  a 
numerical  factor  depending  on  the  lattice  structure  of  the  medium.  A  rough 
estimate  gives  Z  -  10.^  Eq.  (28)  holds  as  long  as  A  is  much  greater  than  the 
lattice  constant  of  the  medium  and  the  spatial  dispersion  of  the  dielectric 
constant  c"  can  be  neglected. ^  Using  Eqs.  (13),  (21),  (22)  and  (28),  we  finally 
obtain  the  total  decay  rates  for  p  in  a  dissipative  dielectric  sphere  as 


^  -  1  +  |  f(w) 
Y0  4 


7  Vy2 >  2 

Re(/e^  £  (2n+l)n(n+l)  » 

n«l  * 
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where  f(u)  *  y^/Yq  is  the  reciprocal  ratio  as  given  in  Eq.  (28).  The  radiative 
rates  can  analogously  be  obtained  from  Eqs.  (26)  and  (27)  by  multiplying  each  of 
them  by  f (u).  The  nonradiative  rates  can  then  be  obtained  again  from  the 
difference  between  the  total  and  the  radiative  rates.  Furthermore,  for  p  located 
at  the  center  of  the  sphere,  we  have  results  analogous  to  those  for  the  case  of  a 

A 

transparent  sphere  obtained  by  Chew  where  for  a  dissipative  sphere  are  given  by 


WWMIUilU'.IUlimiU  ^ 
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(31) 


(32) 


III.  Numerical  Results 

We  have  performed  numerical  studies  for  each  of  the  above  three  classes  of 

problems.  Por  simplicity,  we  treat  one  of  the  two  media  as  a  vacuum  (c  *  1). 

The  other  medium,  if  it  is  transparent  to  the  molecular  emission,  is  taken  to  be 
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glycerol  with  e  ■  2.16,  *  and  if  it  is  dissipative,  it  is  taken  to  be  silver  at 
A  ■  413  run  with  c  *  -4.42  +  0.73  i.*®  Figure  1  shows  the  results  for  a  dipole 

above  a  Ag  spherical  surface  of  radius  a  ■  100  A.  We  see  that  while  the  neglect 

• 

of  the  nonradiative  transfer  (as  in  Chew's  formulation)  leads  to  somewhat  lower 

values  for  y  for  the  radial  dipole  case,  it  leads  to  tremendous  differences  for 

the  tangential  dipole  case  as  r  -*•  a.  This  phenomenon  also  occurs  in  the  flat- 

surface  case,  with  its  origin  due  to  the  parallel  (antiparallel)  orientation  of 

the  "image  dipole"  with  respect  to  the  normal  (tangential)  orientation  of  the 

source  dipole.®  Of  course,  the  extemely  large  values  of  the  rates  at  r  ■  a  are 

not  physical  either,  and  it  is  known  that  more  reasonable  values  for  r  at  r  ■  a 

can  be  obtained  by  considering  the  spatial  dispersion  of  the  dielectric  function 

19 

of  the  substrate  at  close  distances.  Nevertheless,  we  should  add  that  in  all 

the  numerical  calculations  of  Ref.  4,  the  substrate  sphere  is  taken  either  as 

transparent  (c"  ■  0)  or  as  perfectly  conducting  (hence  no  penetration  of  fields 

NR 

and  y  ■  0),  so  that  the  neglect  of  the  nonradiative  retes  does  not  lead  to  any 
error  in  these  cases.* 
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Figure  2  shows  the  results  for  a  dipole  in  a  cavity  of  radius  a  •  4000  A. 
The  outside  medium  is  taken  to  have  a  real  *  2.16.  We  observe  that  both 
radial  and  tangential  dipoles  decay  at  the  same  rate  ("14-  ReE^  (see  Eq.  (23))* 
at  the  center,  with  the  tangential  dipole  showing  appreciable  oscillating 
behavior  as  the  dipole  is  moving  towards  the  cavity  wall.  We  also  note  that  both 
enhanced  and  diminished  values  (with  respect  to  the  free  molecule  case)  for  y  can 
occur  in  this  case.  Figure  3  shows  the  results  for  the  same  cavity  problem  with 
the  outside  medium  being  absorptive  (Ag  at  X  ■  413  nm).  We  notice  that  there  is 
no  oscillating  behavior,  and  the  decay  for  this  case  is  very  small  (unless  the 
dipole  is  far  away  from  the  center)  and  purely  dissipative  in  nature  (we  recall 
that  in  the  exteme  limit  where  the  outside  medium  becomes  perfectly  conducting, 
the  decay  is  zero  as  discussed  in  Sec.  II). 

Figures  4  and  5  show  results  for  a  dipole  inside  an  "aerosol"  of 
transparent  and  dissipative  media,  respectively.  For  the  transparent  case 
(radius  a  ■  4000  A),  we  again  observe  similar  oscillatory  behavior  as  for  the 
flat-surface  case  and  a  possible  "diminution  effect"  for  the  decay  rate.  For  the 
dissipative  case  (radius  a  ■  1000  A),  it  is  interesting  to  note  that  (for  both 
orientations)  most  of  the  decays  are  nonradiative  in  nature,  and  the  surface 
effect  is  hardly  noticeable  except  when  the  molecule  is  very  cloee  to  the 
surface.  This  means  that  most  of  the  molecular  deexcitation  energy  is  dissipated 
within  a  very  small  region  of  the  medium  around  the  molecule,  which  can  be 
understood  from  Eq.  (28)  showing  that  y^  is  very  large  for  molecular  emission 
wavelengths.  To  have  appreciable  surface  eff facts,  one  must  look  for  much  higher 
emission  frequencies,  which  can  occur,  for  example,  in  nuclear  spectroscopy  in 
the  y-ray  range.  Thus  for  our  interest  here  which  focuses  on  molecular 
lifetimes,  the  surface  effect  is  hardly  manifested. 


IV.  Conclusion 


We  have  in  this  paper  presented  a  complete  treatment  of  the  classical  decay 

rates  of  a  molecule  in  the  vicinity  of  a  spherical  surface.  We  have  adopted  a 

different  formulation  of  the  dipole- sphere  interaction  problem  via  the  Hertz 

vector  formalism,  which  is  distinct  from  the  more  common  approach  via  the  Lorenz- 

Hie  expansions  of  the  field  quantities  in  terms  of  complicated  vector  spherical 

harmonics.  We  have  explored  all  the  eight  problems  as  stated  in  the  Introduction 

within  this  classical  approach.  Since  this  approach  is  limited  to  only  two  kinds 

of  mechanisms  (radiative  and  dissipative  transfers)  through  which  the  molecular 

decay  rates  can  be  examined,  some  interesting  features  (e.g.,  possible  enhanced 
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decays  for  a  perfectly-conducting  cavity)  cannot  be  obtained  in  the  present 

formalism.  Recently,  a  quantum  electrodynamical  formalism  (QED)  of  the  spherical 

20 

cavity  problem  was  published,  although  the  electrostatic  image  concept  was 
adopted  and  hence  the  dynamical  nature  of  the  dipole  field  was  omitted.  It  would 
therefore  be  interesting  to  incorporate  the  dynamical  dipole  fields,  as  treated 
in  our  present  paper,  into  the  QED  formalism  to  see  if  it  would  give  rise  to  new 
interesting  phenomena  for  the  decay  rate. 
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Figure  Captions 

1.  Transition  rates  for  the  two  polarizations  of  an  oscillator  molecule  outside 
a  dielectric  sphere  of  radius  a  ■  100  A  for  A  ■  4133  A  and  Cj  *  -4.42  + 
0.731.  The  molecule  is  located  at  a  distance  d  from  the  center,  and  T  and  R 
stand  for  the  total  and  radiative  transition  rates,  respectively.  The  solid 
curve  is  for  the  radial  dipole,  while  the  dashed  curve  is  for  the  tangential 
dipole. 


2.  Transition  rates  for  the  two  polarizations  for  an  oscillator  molecule  in  a 
vacuum  cavity  (Cj  ■  1)  of  radius  a  ■  4000  A  for  A  »  4133  A  and  •  2.16. 

The  molecule  is  located  at  d  from  the  center.  The  solid  curve  is  for  the 
radial  dipole  while  the  dashed  curve  is  for  the  tangential  dipole. 


3.  Same  as  Fig.  2,  except  that  «  -4.42  +  0.731. 


Transition  rates  for  the  two  polarizations  for  an  oscillator  molecule  inside 
a  dielectric  sphere  of  radius  a  ■  4000  A  for  A  «  4133  A  and  ■  2.16.  The 
molecule  is  located  at  d  from  center.  The  solid  curve  is  for  the  radial 
dipole  while  the  dashed  curve  is  for  the  tangential  dipole. 


5.  Same  as  Fig.  4  except  that  tj  ■  -4.42  4-  0.731  and  a  ■  1000  A. 
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